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ABSTRACT 

We study the evolution of the core (r < I kpc) and effective (r < Ug) stellar-mass surface densities 
(El and Eg) in star-forming (SFGs) and quiescent galaxies. As early as z = 3, both populations occupy 
distinct linear relations in log Eg and logEi vs. logM*. These structural relations exhibit almost 
constant slopes and scatter while their normalizations decline with time. For SFGs, the normalization 
in Eg and Ei declines by a factor of ^2 since z = 3. Such mild declines suggest that SFGs build 
dense cores by moving along these relations. We define this evolution as the structural main sequence 
(E-MS), analogous to the star-formation rate main sequence (SFR-MS). Quiescent galaxies follow 
different relations (E^, E^) off the E-MS by having higher densities than SFGs of the same mass and 
redshift. The normalization of E^ declines by a factor of 10 since z = 3, while only a factor of ~2 in 
E)^. Thus, a dense stellar core is present in quiescent galaxies at all redshifts and the formation of such 
core in SFGs is the main requirement for quenching star-formation. Fxpressed in 2D as deviations 
off the SFR-MS and off E)^ at each redshift, the distribution of massive (log(M/M 0 ) > 10.3) galaxies 
forms a universal, L-shaped track that relates two fundamental physical processes: compaction and 
quenching. Gompaction is a process of substantial core-growth in SFGs relative to the evolution in the 
E-MS. This process increases the core-to-total mass and Sersic index, thereby, making “compact” SFGs 
structurally similar to quiescent galaxies. Quenching occurs once compact SFGs reach a maximum 
central density above E^^ = Ei — 0.651og(A'/* — 10.5) > 9.5 M 0 /kpc^. This threshold provides the 
most effective selection criterion to identify the star-forming progenitors of quiescent galaxies at all 
redshifts. 

Subject headings: galaxies: photometry — galaxies: high-redshift 


1 . INTRODUGTION 

Studies of galaxy evolution from the peak of cos¬ 
mic star formation to the present day have matured 
tremendously over the past two decades. The advent 
of large multi-wavelength photometric surveys have en¬ 
abled inferences of the global stellar population prop¬ 
erties such as stellar mass, age, and star formation 
rate (SFR). Large-area surveys such us Sloan Digital 
Sky Survey (SDSS), NMBS, zGOSMOS, UltraVISTA, 
and zFOURGF have robustly established the shape and 
evolution of the mass function of star-forming galaxies 
(SFGs) and quiescent galaxies since z ^ 4, cement¬ 
ing ou r understanding o f galaxy build up an d shut¬ 
down (iPeng et al.l 1201(11: iBra mmer et all 120111 : Ilbert 
et al. 2013: iMuzzin et al.l 1201^ Woo et al. 2013; 
iTomczak et al.l 1201311 . Furthermore, the sensitivity and 
high spatial resolution of Hubble Space Telescope (HST) 
have extended those mass functions further back in cos- 
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mic time (jBouwens et al.ll2010t Finkelstein et al. 2013; 
Oesch et al. 2014), and have made a pivotal contri- 
butioi i to the study of galaxy sizes and morphologies 
(e.g., Ivan der Wei et alTl2012t Shibuya et al. 2015). 
Deep multi-band surve ys with HST, suc h as GOODS 
(iGiayalisco et al.l 1200411 and GANDELS (|Grogin et al.l 
120111 : iKoekemoer et ahl 1201111 have thus provided an 
exquisite dataset to quantify the simultaneous evolution 
of the galaxy stellar populations and structural proper¬ 
ties across cosmic time. 

The consensus is that strong correlations between 
structure and stellar po pulations (i.e., a Hubble se¬ 
quence) exist up to z = 4 (iFranx et al.ll2008l: iKriek et'al] 
l2009al: iWuvts et al.l l2011aL One such relation is the 


star formation rate (SFR) main sec 

uence (SFR-MS) 

(INoeske et al.ll 

20071: 1 

llbaz et al.ll2007l: ISalim et al.ll2007l: 

Pannella et al. 

12009: 

Magdis et al.1 20101: IWuvts et al. 

201 la! Elbaz et all 120111 Rodighiero et al.l 120101: 

Whitaker et al 

.1 120121 iPannella et al.1 I2014H. which is 


thought to describe a relat i vely smooth mode of gal axy 
growth (lElbaz et al.1 120071: iRodighiero et al.l 1201011 in 
which gas i nflow and SFR ha ve reached a steady-state 
phase (e.g., iDekel et'm 120131 1. SFGs on the SFR-MS 
typically have larger sizes and exponential disk profiles, 
while quiescent galaxies of the same mass have more 
concentrated mass profiles (higher Sersic indices) and 
smaller sizes. A dichotomy is also present in the size- 
mass relations, where quiescent galaxies exhibit a much 
steeper slope than SFGs, and a lower normalization, i.e.. 
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highe r densities (iWilliams et al.l l201Ct iNewman et alJ 
120121 : Ivan der Wei et alJl2014li . While both the SFR- 
MS and the size-mass scaling relations evolve with time, 
the fundamental structural differences in SFGs and qui¬ 
escent galaxies are always present, suggesting that hav¬ 
ing concentrated (denser ) surface density profil es is a 
requisite for quenching feauffmann et alJ I2003L 2006 ; 
Schiminovich et all 120071: (Bel. 120081: iCheung et alJl2012l : 
Fang et al.l 120131 : iLang et*^ 1201^ . In other words, 
SFGs must grow dense stellar cores before quenching. 

There is increasing observational evidence that SFGs 
with dense cores exist at every redshift. At z > 2, 
SFGs with the highest central densities are remarkably 
compact and have high Sersic indices and spheroidal 
morp hologies, lacking any signature of an underlyi ng 
disk (IWu vts et a bl 1201 Ibl: iBarro et al l 

Patel et al ■I 120131: Istefanon et al ■I 1201,81: 


an unaerlying 
Ml l2014allR 


Williams et al.l 


201.1 iNelson et al.l 1201-S b These galaxies resemble 


the quiescent population at the same redshift but are 
radically different from other SFGs that have irreg¬ 
ular and clumpy a ppearances (lElmegreen et al.l 12004 
iGenzel et al.l 120081 : iGuo et all I2015D . This suggests 
that compact SFGs are formed by strongly dissipa- 
tional processes. Some of these processes, like merg¬ 
ers and disk instabilities, are indeed expected to be 
more frequent at earlier tim es due to the higher gas- 
to-total mass ratio s in SFGs (|Tacconi et al.l[20ll . 120131 : 
iDaddi et al.l 12010( 1. The increased gas mass relative 
to the SFR makes such systems prone to gravitational 
collapse on scales of ^1 kpc, causing substantial core¬ 
growth resulting from a gas-fed cent ral starburst and/or 
an inward migratio n of clumps (iDekel et al.l l2009bl: 
Geverino et al.l 120101: Genel et "Si 12014 : Wellons et al l 


2014 iGeverino et al.! 120151 : iZolotov et al.l 1201,^ At 

lower redshifts, SFGs with dense stellar cores have 
clearly recognizable disk structures, but their profiles 
are dominated by a central bulge (IWuvts et al. [2M^ 


iBruce et al.ll20l3 iBruce et al.ll2014ai : iLang et al.l ^14) . 

Interestingly, quiescent galaxies at low-z also seem to 
have bulge-bdisk morphologies (iM cGrath et al.l 120081 : 
iBundv et al.ll2010l : | van der Wei et al .1120111 ) . suggesting 
that quenching takes place among galaxies with similar 
morphologies. 

The common denominator in the evolution of massive 
galaxies described above is the growth of a dense stel¬ 
lar core. This suggests that it is possible to describe the 
general processes of structural growth and star-formation 
quenching using a unique quantity tracing central stellar 
mass density. Here we study the evolution since z ~ 3 
of the stellar mass surface density within a radius of 
1 kpc, El. We build upon previous results at lower red¬ 
shift which show that Ei is closely related with quies¬ 
cence, and follows a much tighter correlation with stellar 
mass than the effective radius or the effective surface 


density, Ilpfe.g., ICheunff et alJ 12012 

. iFanff et al.l 120131: 

Ivan Dokkum et al.l l2014l iWoo et al.l 

2015[). We aim to 


answer whether this relation holds at high redshift, if it 
is a more fundamental quenching predictor, and if the 
global build up and quenching of SFGs can be described 
in simple terms using Ei, i.e., if it can be used to track 
galaxies in transit from the star-forming to quiescence 
phase. 

The paper is structured in three parts: In § [2l we give 
an overview of the observational dataset and the esti¬ 


mated properties of the galaxy sample. In § 13.11 we ana¬ 
lyze the correlations in effective and central mass sur¬ 
face densities, Ee and Ei, vs. stellar mass for SFGs 
and quiescent galaxies from z = 3 — 0.5. In§ 13.21 we 
study how the best-fit relations in Eg and Ei relate to 
each other for galaxies with a Sersic mass profile. In 
§ 13.41 we analyze galaxy evolutionary paths based on the 
best-fit Eep structural relations and predictions of hydro- 
dynamical simulations. We find that the structural evo¬ 
lution of SFGs can be described in terms of 2 phases: 
A steady phase of size and core growth that we call the 
E “main sequence ”, an d a “compaction” phase of rapid 
core growth. In § 13.51 we study the relative distance of 
massive SFGs and quiescent galaxies from the SFR-MS 
and the quiescent structural relations as a function of 
redshift. We find that their distribution defines a univer¬ 
sal (redshift-independent) L-shape sequence that traces 
2 fundamental transformations: the compaction of SFGs 
to form a dense stellar core, and the quenching of star- 
formation in those galaxies at maximum central density. 
Lastly, in S I3.6l we discuss the evolution in the structural 
properties and visual morphologies of galaxies within the 
compaction-quenching sequence as a function of time. 

Throughout this pa per, we quote magnitudes in the 
AB system, assume a (|Ghabrierll2003l) initial mass func¬ 
tion (IMF) and adopt the following cosmological param¬ 
eters: h) = (0.3, 0.7, 0.7). 

2. DATA 

This paper is based on a sample of massive galaxies 
built from the HST/WFCS F16 0W selected cata log for 
the GANDELS GOODS-S field (IGuo et a,1 . 112011 ) . Gon- 
sistent, m ulti-wavelength pho tometry was measured us¬ 
ing TFdTJ^aMle^tlayllOOl) , implemented as described 
by IGuo et al.l (|2012DTPhotometric redshif ts were com- 
puted using EAZY ((Brammer et al.l 120081) and yielded 
errors of Az/(1 -|- z) = 3% . Stellar masses were de¬ 
rived using FAST (iKriek et al.ll2009b[) and based on a 
grid oflBriizual fc Gharlotf i 20031 ) models that assume a 
iGhabried (I2003D IMF, solar metallicity, ex ponentially de¬ 
clining star formation histories , and the lGa,lzetti et al.l 
(|2000i) dust extinction law (see ISantini et al.ll2015l for a 
detailed description). 

Star formation rates (SFRs) were computed by com¬ 
bining IR and res t-frame UV (un corr ected for extinc- 
tion) luminositie s (iKennicuttI 119981 and iBell et ani2005D 
and adopting a iGhabrieil (|2003l) IMF (see iBarro et al.l 
120111 for more details): S'Fi?uv-i-iR = 1-09 ^ 10 “^°(Lir-|- 
3.3L28oo)- Total IR luminosities ( Atr= L[8-1000/rm]) 
were derived from iGharv fc Elbad (|2001l) templates fit¬ 
ting MIPS 24Mm fluxes, app lying a Herschel-hased re¬ 
calibration (lElbaz et al.ll201ll ). For galaxies undetected 
by MIPS below a 2cr level (20p,Jy) SFRs come from rest- 
frame UV luminosities th at are corrected for extinction 
as derived from SED fits (|Wuvts et al.ll20iral) . 

The half-light radii, measured along the major axis, 
and Sersic indices were determined from HST/WFC3 H 
images using GALFIT (|Peng et al.ll2002l ) with PSFs cre¬ 
ated and pro cessed to replicate the co nditions of the ob¬ 
served data ((van der Wei et ai]l2012D . The stellar mass 
profiles were computed by fitting multi-band SEDs de¬ 
rived from surface brightness p rofiles in 9 HST bands 
measured with IRAF/ellipse (see lLiu et al.l[2M^ and Liu 
et al. in prep for more details). Following IWuvts et al.l 
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Fig. 1. — Star-formation rate vs. stellar mass for galaxies in the CANDELS/GOODS-S field at 0.5 < z < 3.0. The blue lines show 
our fits to the S FR-M S at each redshift, which agree with previous results above log(M/M0)> 10. The gray lines show the results of 
IWhitaker et al.l 1)20131) . which highlight the change in the slope of the SFR-MS at lower masses. SFGs are indicated by the blue points. 
Quiescent galaxies (red points) are selected to lie 0.7 dex below the SFR-MS, i.e., below the red dashed lines. 


(IMl , we impose an additional constraint on the spa¬ 
tially resolved SED-fit by requiring that the integrated 
profile matches the observed flux in IRAC chi and ch2. 
We apply this constraint by assuming that the integrated 
IRAC-F160W color is the same at all radii. 

The multi-band HST mosaics were PSF-matched to 
the resolution of F160W, that has a half width at half 
maximum of HWHM=0.”09. The profiles have an in¬ 
trinsic spatial resolution ranging from r — 0.6 — 0.7 kpc 
within the redshift range of the sample and thus resolve 
the inner 1 kpc of the galaxies. However, part of the 
light can be smeared to larger radius. To correct for 
this effect, we estimate a Sersic dependent correction to 
the mass profile within 1 kpc. The correction is com¬ 
puted from a grid Sersic profiles with n = 0 — 4 and 

re = 0.1-10 kpc degraded to F160W resolution at 

different redshift. The correction ranges from ^ 0.4 dex 
at n > 2 to ^ 0.2 — 0 dex at n = 1 — 0. 

Our goal is to analyze the distribution of massive SFGs 
and quiescent galaxies at 0.5 < 0 < 3 in the structural 
scaling relations in order to understand their differen¬ 
tial evolution. The deep CANDELS photometry enables 
the selection of stellar mass complete samples down to 
log(M/MQ)^ 10.3 at 2 : = 3 (e.g., Tal et al. 2014). 
Thus, we use this threshold to limit our sample selec¬ 
tion for the analysis in § 3.4 and 3.5. Nonetheless, in 
the next section, we use a larger sample of galaxies down 
to log(M/MQ) = 9 (~70% complete) to characterize the 
SFR and structural correlations over a wider dynamical 
range. 


3. RESULTS 

3.1. SFR and structural scaling relations 

Figure [T] show the distribution of SFR vs. stel¬ 
lar mass for the galaxies in GOODS-S. The major¬ 
ity of star-forming galaxies follow a relatively tight re¬ 
lation between SFR and stellar mass. The observed 
relation suggests that galaxy star formation histories 
are predominantly regular and smooth, i.e., galaxies 
grow in a secular mode which is us ually referred to as 


the SFR main sequence (SFR-MS: INoeske et al.l 

2007; 

Elbaz et alJl2007l:ISalim et alJl2007l: 

Pannella et al.l 

2009; 

Magdis et al. 20101: Wuvts et al.l 

201 lat Elbaz et al. 

20111: IRodighiero et al.l 120101: IWhitaker et al.l 120121: 


iPannella et al.l 120141) . Following previous works, we 
characterize the SFR-MS as a single power-law at 
log(M/M0) >10. 


log SFR = /r 



10.5 


+ logC 


( 1 ) 


Note that at l ower masses, the S FR-M S e xhibits a steeper 
slope (e.g., IWhitaker et al.l (|2014ll : iSchreiber et al.l 
(j201,5l lb The best-fit SFR-MS at ever y redshift is in 
excelle nt a greement with the res ults of IWhitaker et aH 
(|2014f) and lSchreiber et ^ (120151) at the high-mass end. 
The slope and the normalization are consistent with their 
values within the errors (Table[T|). We select star-forming 
galaxies above and quiescent galaxies below a threshold 
of ASFR = —0.7 dex (dashed red line), where ASFR = 
log SFR — log SFR'^®. This classification line is ^ 2tT be¬ 
low the SFR- MS, which has a typical observational scat¬ 
ter of 0.3 dex (IWhitaker et al.ll201-S iSoeagle et al.ll20l4 
ISchreiber et al.ll2015l) . 

Panel A of Figure H] shows the distribution in effec¬ 
tive surface mass density, Eg = 0.5M*/7rrg, vs. mass for 
the galaxies in Figure [TJ SFGs and quiescent galaxies 
follow well-defined size-mass relations, which are char¬ 
acterized byMowdinea£ 2 £l 3 iioSSjJ,OK££_o^.alogJW(e^, 
Law et al.ll201^lMosleh et al.ll2012l:ISzomoru et al.ll20i2l : 

Newman et al.ll2012t Ivan der Wei et al.ll2014l) . Those re¬ 

lations can be expressed in terms of Eg: 



logA( 2 :) 


( 2 ) 


where a is related to the slope of the size-mass rela¬ 
tion a as a = 1 — 2a. The red and blue lines show 
the best-fit Eg relations for quiescent and star-forming 
galaxies. The best-fit slopes are relatively constant with 
time, ^ —0.5, ^ 0 . 6 and agree well with the 

results of Ivan der Wei et al.l (j2014[ ) for the size-mass re¬ 
lations of both populations (a®^ ^ 0.2 and a^ ~ 0.8). 
There are too few quiescent galaxies at 2.2 < z < 3.0 
to accurately fit the slope, and so we fix the slope (but 
not the normalization) in this bin to match the value at 
1.4 < z < 2.2. Note that quiescent galaxies have a steep 
slope in the size-mass relation, which leads to an anti¬ 
correlation in E^. Meanwhile, SFGs have a shallow slope 
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Fig. 2.— Surface density vs. stellar mass for galaxies in the CANDELS/GOODS-S field at 0.5 < 2: < 3.0. Panels A and B show the 
surface density within the effective radius, Eg, and within the inner 1 kpc, Ei, respectively. The blue and red circles show SFGs and 
quiescent galaxies selected using the SFR criterion of Figure^ The thick blue and red lines depict the best-fit log E — log M* relations for 
the two populations. The dashed red lines show the 2a scatter around the quiescent relation. SFGs and quiescent galaxies exhibit clear 
and distinct scaling relations since 2 : 3, which are well-described by single power-laws. The scatter in the Ei relations is a factor of 2 

smaller than that of Eg. The slopes of the scaling relations remains approximately constant with time. The zero-points of the star-forming 
relations decline slowly with time (see Figure[3jl. Based on this smooth evolution, we speculate that SFGs follow, on average, evolutionary 
paths along that relation. We define this track as a structural “main sequence” (E-MS) . At log(M/MQ)> 10 we find an increasing number 
of compact SFGs with high surface densities, similar to those of quiescent galaxies (green circles selected within the quiescent relation in 
El). We capture this trend by fitting a second-order polynomial to logE — logM*, which shows a steeper slope at the high-mass end (thin 
blue l ine). T hese galaxies may deviate upwards from the E-MS due to dissipational “compaction” events that cause a rapid core growth 
(see ^ 13. 4. 211 . Fit parameters are given in Tables^ and [2 


in size-mass (a®^ < 0.5) and a positive correlation be¬ 
tween Eg’' and mass. The scatter in Eg is consistent with 
^2x that of the size-mass relations, cr(logEe) ~ 0.5 dex 
and 0.3 dex, for SFGs and quiescent galaxies, as expected 
from AlogEe oc 2Alogre. The redshift-dependent nor¬ 
malizations decline from z = 3 to z = 0.5. Such decline is 
much steeper for quiescent galaxies than for SFGs (1 dex 
vs. 0 . 3 dex) as noted in previous works (iBuitrago et al.l 
I2008t INewman et al.l[20T^ Ivan der Wei et al.ll2014ll . 

The bottom row (B panels) of Figure [5] shows the red- 
shift evolution of the central surface mass density within 
1 kpc, Eikpc = M{< 1 kpc)/7r(l kpc)^, versus the stel¬ 
lar mass. Similarly to Eg, we characterize the observed 
correlation in Ei as a log-linear relation: 



10.5 


-blogB(z) 


( 3 ) 


Again, we find clear correlations for both SFGs and qui¬ 
escent galaxies at every redshift since z ^ 3. The slopes 
of these relations are positive and relatively constant 
with time, /3®^ = 0.9 and /3^ = 0.7. By comparison with 
Eg, the dispersion is ~ 2x tighter, cr(logEi) ~ 0.25 dex 
and 0.14 dex, for SFGs and quiescent galaxies , in good 
agreement with the results of iFang et al.l (I2013D at z = 0. 
The normalization of the star-forming Ef’' relation de¬ 


clines by ~0.3 dex from z = 3toz = 0.5 similar to the 
evolution in E®’'. Interestingly, for quiescent galaxies, 
E^ declines by a similar amount, in stark contrast with 
the strong decline of ~ 1 dex in E^. The lower scat¬ 
ter and weaker redshift evolution indicates that Ei is a 
more robust and reliable structural parameter than Eg. 
Note also that, by measuring the mass inside a fixed 
physical core aperture, Ei is closer to the concept of a 
cosmic clock, i.e., it only increases with a stellar mass 
growth, unlike the effective size (re), which can also de¬ 
crease (e.g., due to a substantial mass growth closer to 
galaxy center, or to fading of a extended star-forming 
region). 

Two main conclusions arise from the distribution of 
SFGs and quiescent galaxies in Figure [2] 1) the slopes 
of the structural relations are almost constant with time, 
and the normalizations decline for all of them; 2) at any 
redshift, quiescent galaxies are denser than SFGs of the 
same mass, although this difference declines with time 
in Eg due to fast evolution of E^. A reasonable assump¬ 
tion based on the first conclusion is that the evolutionary 
paths of individual SFGs approximately follow the best- 
fit Eg’^relation which, by analogy with the SFR-MS, de¬ 
fines a structural main sequence, E-MS, that could be 
interpreted as phase of smooth structural growth. The 
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TABLE 1 
Power Law Fits 



SFRms- SFGs 

Ee — Quiescent 

Se - SFGs 

El — Quiescent 

El - SFGs 

redshift range 

P logC 

a logA 

a logA 

/3 logB 

d logB 

0.5 < 2 < 1.0 
1.0 < 2 < 1.4 
1.4 < 2 < 2.2 
2.2 < 2 < 3.0 

0.19 ±0.08 1.21 ±0.02 
0.53 ±0.07 1.44 ±0.04 
0.64 ±0.06 1.75 ±0.05 
0.68 ±0.06 1.92 ±0.05 

-0.42 ±0.13 9.15 ±0.06 
-0.45 ±0.14 9.53 ±0.05 
-0.52 ±0.14 9.91 ±0.07 
-0.48 ±0.17 10.16 ±0.05 

0.60 ±0.04 8.46 ±0.04 
0.60 ±0.05 8.54 ±0.05 
0.64 ±0.05 8.68 ±0.04 
0.56 ±0.06 8.80 ±0.04 

0.65 ±0.03 9.53 ±0.05 
0.65 ±0.04 9.64 ±0.04 
0.64 ± 0.03 9.74 ± 0.05 
0.67 ±0.04 9.80 ±0.05 

0.89 ± 0.03 9.12 ±0.03 
0.88 ± 0.03 9.16 ±0.04 
0.86 ±0.04 9.25 ±0.03 
0.89 ±0.04 9.33 ±0.05 


Notes, power law coefficients parameterizing the evolution of the SFR (logSFR= /ilog(M* — 10.5) + logC), the effective 
mass surface density (logEe= Q:log(M* — 10.5) + logA) and the mass surface density within the central 1 kpc (logEi = 
/31og(M* — 10.5) + logB) relations with stellar mass as a function of redshift. 


TABLE 2 

2ND-ORDER Power Law Fits 



Se-SFGs 

El-SFGs 

redshift range 

Oil a2 logA 

di d 2 logB 

0.5 < 2 < 1.0 
1.0 < 2 < 1.4 
1.4 < 2 < 2.2 
2.2 < 2 < 3.0 

0.09 ± 0.08 1.10 ± 0.07 9.25 ± 0.03 
0.15 ±0.07 1.19 ±0.06 9.25 ±0.02 
0.16 ±0.06 1.18 ±0.07 9.36 ±0.02 
0.14 ±0.06 1.10 ±0.05 9.41 ±0.02 

0.08 ± 0.04 0.80 ± 0.03 8.56 ± 0.05 
0.08 ± 0.03 0.95 ± 0.04 8.65 ± 0.04 
0.09 ±0.03 1.12 ±0.03 8.82 ± 0.05 
0.11 ± 0.04 1.10 ± 0.04 8.93 ± 0.05 


Notes. Power law coefficients parameterizing the evolntion of Ee and Ei vs. 
stellar mass as 2nd order polynomials: logEe= ai (log M)^ + Q 2 log M + log A 
and logEi= /3i(logM)^ + P 2 logM + logB 


second conclusion indicates that increasing the surface 
density by forming a dense stellar core is a requisite for 
quenching the star-formation, as similarly suggested in 
previous works. Together, these conclusions suggest the 
need for an intermediate phase to bridge the E-MS and 
quiescent sequence, in which SFGs become more com¬ 
pact and centrally concentrated before quenching. We 
further discuss the possible evolutionary paths of SFGs 
and quiescent galaxies in § 13.41 


3.2. Surface mass profiles: relating Eg and Ei 

In this section we analyze the relation between the 
Ee and Ei structural relations for SFGs and quiescent 
galaxies in terms of their surface density profiles. Galax¬ 
ies’s surface bright ness and stell ar mass profiles are typ¬ 
ically described bv iSersid (I1963I1 profiles, in which case, 
the central and effective densities are related to one an¬ 
other by the profile parameters. Assuming the following 
characterization of the mass profile: 


M{r) 


M, 


exp 


( - b„ 



( 4 ) 


where n is the sersic index and Mg is the effective mass 
at Te, El and Ee can be determined from one of them 
for a given value of n and r^. Note that in doing so we 
assume that the mass and light profiles are the same, 
Te.mass = Te,light- This is a better approximation for qui- 
escen t galaxies, which exhibit only weak color gradients 
(e.g.,[^o et al.l[2012l: iSzomoru et al.l[2M^ iWuvts et al.l 

I 2 OI 2 D . Nonetheless, the parametrization is useful to dis¬ 
cuss the expected differences in Ei and Ee when the as¬ 
sumption fails. 

By integrating Equation 3] we obtain: 


log Eikpc - log M* = - log TT - log 7(2n, h^r^ ^Z"-) (5) 


where 7 is the incomplete g amma function, whic h de¬ 
pends on n and Vf. (see e.g., iGraham et al.l 120051) . As 


shown in Figure SI the 7 function can be approximated 
as a second order polynomial for a given n. 

log 7 ( 2 n, 6„r“^/") = cq -k Ci log re -k C 2 (log re)^ (6) 

where the Ci coefficients depend on the Sersic index. 
Combining the linear term of this equation with the 
best-fit relations for Ei and Ee as a function of mass, 
and the definition of the latter in terms of re (log re = 
0.5(log M*—log Ee)), we obtain the following relations for 
the slope and the scatter of the logEi — logM* relation 
in terms of the slope in the Ee and size-mass relations. 

/3~ 1 -k 0.5ci(a — 1); cr(log Ei) ^ 0.5cicr(log Ee) (7) 

/3~l-cia; tT(logEl) ^ Cl cr(log re) 

Quiescent gal axies typically have hi gh Sersic values , 
n = 3 — 4 (iWilliarns et al.l 120101: Bell et ahl 1201211 . 
Therefore, for a value of ci(4) = 0.55, we obtain 
tT(logEi) ^ 0.3(T(logEe), which explains the smaller 
scatter in Ei with respect to Eg. Furthermore, for a 
value of a ~ —0.5, the expected slope of the Ej* relation 
is /3 ^ 0.6, which is consistent with the best-fit value. 
Thus, a single Sersic model is good approximation to ex¬ 
plain the differences in the E^ and E^ relations. 

SFGs exhibit a broad range in Sersic values, and their 
profiles can often deviate from a single Sersic profile , 
e.g., in a disk wi th a denser bulge (|Wuvts et al.ll20llal : 
iBruce et al.ll2012D . The color code in Figured illustrates 
the spread in Sersic values for SFGs in the Ei vs. mass di¬ 
agram. At a given stellar mass, the Sersic increases with 
El. This correlation suggest that as SFGs grow a dense 
stellar core, their mass profiles become more centrally 
concentrated. I.e., the width of the Ei distribution for 
SFGs depends on their structural properties. The SFGs 
in the overlapping region with Ej, in particular, exhibit 
Sersic values similar to those of quiescent galaxies. The 
blue and cyan lines show the predicted Ef*’ relation based 
on Eg'’ for n = 1 and n = 2. Qualitatively, the predic- 
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redshift 

Fig. 3.— Redshift evolution in the normalization of the Si (stars; 
logB) and Se (circles; log A) relations for SFGs (blue) and quies¬ 
cent galaxies (red). The solid and dashed lines show the best-fit 
relations to the normalizations in Si and Se as a function of red¬ 
shift. The purple star shows the zero-poi nt of the quiescent Si re¬ 
lation at z = 0 from I Fang et al.l 1120131 1 . The dashed gray line 
shows the predicted evolution in Sj due to minor mergers, assum¬ 
ing that Sj evolves only due to mass growth outside 1 kpc. The 
dashed black line shows the predicted evolution in Sj assuming 
that S^ evolves only due to mass loss and adiabatic expansion. 
The gray circles show the predicted normalization in S^ inferred 
from S^ using Equation [5] and a S&sic of n = 4. 

tion for n = 2 matches the slope and normalization of 
the best-fit S®'' relation (gray). However, the relations 
determined from E®*' with n < 2 have lower normaliza¬ 
tions than the observed distribution of SFGs with similar 
Sersic values. This suggests that the normalization of the 
Eg distribution for SFGs also increases with n. 

3.3. Redshift evolution of the normalization of the 
quiescent structural relations 

Using equation [5] and the relation between Eg and rg, 
we quantify the evolution in the normalization of the 
structural relations (log A and log B) due to different evo¬ 
lutionary processes that preserve the Sersic profile but 
cause a relative change in size and stellar mass character¬ 
ized as A log Ve = rjA log M. If the slope of the structural 
relations remain constant in the process, the normaliza¬ 
tions change as AlogEg^i — (a,/3)A logM, which leads 
to: 

Alog A = (1 — 2?? — a)AlogM (8) 

A log B (1 — ciT] — P)A log M 

These relations are useful to test whether some of the 
proposed evolutionary processes are consistent with the 
simultaneous evolution in the normalizations of both 
structural relations (Table [T|). For quiescent galax¬ 
ies, where the relative evolution in the normalization 
of Eg vs. El is largest, there are three main evolu¬ 
tionary channels that can explain such large difference. 
The leading explanation is that quiescent galaxies ex¬ 
perience a large increase in size compared to the mass 



Fig. 4. — Ratio of the mass in the central 1 kpc to the total mass 
(defined as C/T) as a function of the effective radius, re, and n for a 
single-Sersic mass profile. The C/T ratio follows a gamma function 
(Equation O, which can be expressed as a second-order polynomial 
in Ve for different values of n (bottom-left corner; Equation in. 
The linear coefficient of those relations, ci, determines the relative 
values of the slope and typical scatter of the Eg, Si and Ve scaling 
relations with stellar mass (Equation O . The best-fit structural 
relations for quiescent galaxies, EJ and agree well with one 
another for a Sersic profile of n = 4 (see gray and red lines in 
Figure [Sjl . 


growth due to mi nor mergers (i.e. , accr e tion of smaller 
satell i te galaxies : iBezanson et abl 120091 : iHookins et al.l 
I2009I lOser et iall 120121 ). Alternatively quiescent galax¬ 
ies can have a large increase in size {puff up) due to 
feedback or stellar mass loss associated with passive evo¬ 
lution (i.e., death of old stars) which cause adiabatic 
expansion (iDamianov et aI1l20lH iPoggianti et al.l[2013l : 

IPorter et al.ll2014Hvan Dokkum et al.ll2014l) . Lastly, the 

strong decline in E^ could be caused by the arrival of 
new quiescent with progressively larger sizes at lower 
redshifts fe.g.. IPoggianti h al.ll2013l iCarollo et al.ll20ll : 

iPorter et al.l 12014 ). Each of these scenarios leads to a 
different evolution in the normalization of E^ and E)^. 
Attending to the best-ht relations in the upper-left cor¬ 
ner of Figure [31 the expected evolution in the ratio of 
normalizations is A log A/A log B ^ 5. 

If the evolution of quiescent galaxies is driven by mi¬ 
nor mergers, the expectation is that the core mass re¬ 
mains relatively unchanged in the process. Thus the 
normalization in E)^ changes only due to the total mass 
growth as A logB ^ —0.65A log A/, while the normaliza¬ 
tion in E^ follows Equation 8 as A log A ^ —l.TAlogM, 
for c onservative value of 77 ^ 1.6 (e.g., iNewman et al.l 
I 2 OI 2 D . The predicted relative change in such case is 
AlogA/AlogB ^ 2.6, which is substantially lower than 
the observed ratio (gray line in EigurejS]). Therefore a 
simple minor merger scenario is inconsistent with the ob¬ 
served results. This tension can be partially alleviated if 
mergers cause a larger size growth (77 ^ 1.6 ) or if they 
also cause an increase in Ei, e.g., due to projection ef¬ 
fects in the 2D surface density . The latter, however, is in - 
consistent with the results of Ivan Dokkum et al.l (I2014D . 
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Fig. 5.— Si vs. stellar mass for only SFGs at 1.4 < z < 3.0 
color coded by median Sersic index. The gray lines show the best- 
fit log Si — logM* relations for SFGs and quiescent galaxies at 
z = 2 depicted in Figure [2] The SFGs in the overlapping region 
with the quiescent Si relation (dashed gray lines) exhibit higher 
Sersic index. This could indicate that those SFGs arrive on the 
quiescent sequence as a result of a structural transformation that 
makes their mass profiles more centrally concentrated. The red 
line shows the Sj® relation inferred from the observed S^ relation 
for quiescent galaxies using Equations [5] and |6]and Sersic value of 
n = 4. The observed and predicted Si relations show an excellent 
agreement. The blue and cyan lines show the predicted S®^ re¬ 
lation for SFGs derived from S®®^ assuming n = 1 and n = 2, 
respectively. Qualitatively, the slopes agree with the observed re¬ 
lation, but the systematically lower normalizations suggest that 
the zero-point of S®'^ depends on the Sersic. 

who showed that the number density of the most massive 
cores decreases with time. 

An alternative, or most likely complementary, scenario 
to explain the fast decline in the normalization of is 
that newly quenched galaxies at lowe r redshifts have 
large r sizes and lower den sity cores le.g.. lPoggianti et'm 
120131 : l^rollo et ^120131) . Assuming that all those qui¬ 
escent galaxies have also high Sersic indices (n 4), we 
use equations 5 and 6 estimate the redshift dependent 
normalization in based on the observed values of the 
the normalization in (gray circles in Figure [3]). This 
prediction exhibits a much better agreement with obser¬ 
vations down io z = 0.75, suggesting that a single Sersic 
model is a good approximation at those redshifts. The 
main difference with a minor merger scenario is that the 
latter increases the size and preserves the core density 
at the expense of breaking the single Sersic mass pro¬ 
file. Note however that at z = 0, the predicted value of 
based on the E^ value of (|Fang et al.ll2013ll deviates 
from the expected trend towards larger sizes for a given 
central density. This could signal a more prominent role 
of minor mergers at lower red shifts or a change in the ra¬ 
tio of wet-to-dry m ergers te.g. lLopez-Saninan et al.ll2012l : 
iPorter et al.ll201^ . 

Another evolutionary scenario that preserves the shape 
(Sersic index) of the mass profile is adiabatic expansion. 
In this case, galaxies puff up due to a decline in the 
gravitational potential caused by the death old stars. 


For a value of ry = —1 (e.g., iDamianov et al.l [200911 . 
E^decreases as A log A ~ 3.5AlogM*, while E^ de¬ 
creases as AlogB ^ AlogM*. Thus, the ratio between 
the two is A log A/AlogB ~ 3.5, which is larger than 
the prediction from minor mergers, but still underesti¬ 
mates the observed trend (dashed black line). The nearly 
unity relation in log B with stellar mass indicates that the 
mass loss required to reproduce the observed evolution 
Q in E)^ from adiabatic expansion is A l og AA (1 -|- 
(Figur e [3|). However, as noted by Ivan Dokkum et alJ 
()2014f l. a simple formation model in which quiescent 
galaxies quenched as early as z 5 implies a much slower 
mass decline (^ (1 -I- z)° °®). 

As pointed out in previous works, the most likely sce¬ 
nario is that all the evolutionary channels above play a 
role in the evolution of the normalizations. Nonetheless, 
quantifying their relative contribution as a function of 
time requires precise estimates of th e ages and quench - 
ing times of quiescent galaxies fe.g.. iBelli et al.l[2014bll . 
as well as large enough samples to characterize the num¬ 
ber densities at the extremes of the distributions (i.e., the 
smallest or densest galaxies) w hose disappearance indi- 
cates the nee d for a size growthlyan der Wei et al.ll201^ 
or mass loss ((van Dokkum et al.ll2014r i~ 

In passing, we note that the decreasing trend with 
tim e in Ei has also im plications for the galaxy dynam¬ 
ics. iFang et all (|20I3fl found an excellent correlation in 
El and the velocity dispersion in the central I kpc, which 
closely followed the virial theorem, Ei oc cr^. If this cor¬ 
relation holds at high-z, the decrease in Alogi3(z) ^ 
0.3 dex from z = 3 to z = 0 implies that the Faber- 
Jackson (jFaber fc JacksonllI976ll relation at z > 2 should 
be 0.15 dex higher than the local value, as suggested in 
IBelli et al.l (l2014lbl) . 

3.4. Galaxy evolutionary paths in E, Ei vs. M*; the 
structural “main sequence” and compaction: 

In this section we study possible evolutionary paths for 
SFGs and quiescent galaxies in Ee and Ei vs. mass based 
on the evolution in the normalization of the structural 
relations computed in S 13. H and predictions of theoretical 
simulations. 

3.4.1. Observational trends 

Based on the conclusions in § 13.11 we speculate that 
the relatively constant slope and weak evolution of the 
normalization in E®'’ and Ei suggest that SFGs follow 
evolutionary paths along these sequences, which we de¬ 
fine as the E-MS. In the following discussion we adopt 
this assumption, and thus we refer to the E-MS as log- 
linear tracks in AlogE(ei) = {a, l3)Alog M with slopes 
of a ^ 0.7 and /3 ^ 0.9. 

SFGs in the E-MS increase their surface density with 
time. However, since the slope of the E®*' relation is 
,5 < 1 the mass in the core grows more slowly than the 
total mass of galaxy, i.e., the core-to-total mass ratio de¬ 
creases with time. This is consistent with the notion of 
inside-out growth of an exponential profile (disk), due 
to galactic-scale star-formation and accretion of higher 
angular momentum material, which causes both Vg and 
El t o increase proportionally to the stellar mass growth 
(e.g.. lNelson et al.ll20I2ll20I3l) . This scenario also agrees 
with Figure [5l which shows that the E-MS (gray) de¬ 
scribes an evolution at approximately constant Sersic. 















































Fig. 6.— Galaxy evolu tionary tracks in Se (le ft) a nd Si (rifilit) v s . stell ar mass as a function of time for the Vela set of 35 hydrodynamical 
simulations described in ICeverino et J] II2014I1 and IZolotov et aP II2015I) . The black and gray lines show the median evolutionary tracks 
for galaxies above and below a threshold of log(M/M0)= 10.2 at 2 = 2. The dashed and solid parts of the tracks indicate the evolution at 
redshift higher and lower than 2 = 2.5, respectively. The blue and red contours show the distribution of SFGs and quiescent galaxies at 
2 = 1 — 2. The blue and red lines show the best-fit Se and Si relations at 2 ~ 1 and 2 ~ 2 from Figure[2] The model tracks are in excellent 
agreement with the observed distributions in Se and Si (see also Figure 12 of Ceverino et al. 2015). Moreover, the tracks exhibit two 
fundamental phases with different slopes (blue arrows): (1) the S-MS, a phase of smooth structural growth that follows the best-fit Se,i 
relations from the data (a ~ 0.6, 0 ~ 0.9), and (2) compaction, a period of steep core-growth (a ~ 1.5, 0 ~ 1-3), usually triggered by a 
strongly dissipational event. The simulations exhibit a downsizing trend, such that the most massive galaxies evolve earlier and experience 
a stronger compaction event (open star) due to higher gas fraction at high-z. The different loci of SFGs and quiescent galaxies is caused 
by the downsizing trend: massive galaxies form the backbone of the sequence by quenching first (and rapidly) at high-z, while low mass 
galaxies arrive on the sequence from below in a late compaction event. The green line shows the track of Velal3, which exhibits the latest 
compaction at 2 ~ 1, remaining the longest in the S-MS . The red arrows indicate possible evolutionary tracks for quiescent galaxies for a 
minor merger or adiabatic expansion scenario (see Ij 13.311 . 


For quiescent galaxies, which lack in-situ star- 
formation, the evolutionary paths are thought to be 
driven by either minor mergers or adiabatic expansion 
due to mass loss (see previous section). In a minor merger 
scenario, quiescent galaxies follow flat tracks in Si and 
steeply declining tracks in Se (a ^ —2.2). In an adia¬ 
batic expansion scenario, quiescent galaxies lose stellar 
mass and puff up, both within the Vf, and the central 
I kpc, which causes a steep decline in Se (a = —3) 
and Si (/3 = —I). The evolutionary tracks for quies¬ 
cent galaxies in these scenarios are shown as red arrows 
in Figure m 

Next, we focus on the evolutionary paths that take 
SFGs to the higher-density structural relations of qui¬ 
escent galaxies. To first order, SFGs growing along 
the S-MS intersect the quiescent structural relations at 
log(M/MQ) > II.5. However, confining quenching to 
this evolutionary track would overp roduce the number 
ofmassive quiescent galaxies (e.g., Ivan Dokkum et aP 
This unique track also requires a population 
of extremely dense, low mass SFGs that are not ob¬ 
served. An alternative scenario to explain the emer¬ 
gence of quiescent galaxies is that some SFGs follow a 
steeper path upwards from the E-MS as result of some 
period(s) of fast core growth. These periods can be 
caused by strongly dissipational, “comp action” events, 
e.g., m ajor mergers (iH opkins et al.ll200 8D or disk insta¬ 
bilities (jElmeereen et al. 120081 : iDekel et al1l2009b[ ) which 
funnel large amounts of gas to the center of the galaxy. 
Gompaction events enhance star-formation in the core 


increasing the central density and, in the most extreme 
cases, collapsing the whole galaxy to a much smaller ra¬ 
dius. 

The c omp action scenar i o was discussed in iWnvts et al.l 
(|20IIbD and iBarro et ^ (|20I3D on the basis that some 
massive SFGs seem to lie on the steep size-mass rela¬ 
tion of quiescent galaxies. These compact SFGs exhibit 
smaller sizes and higher Sersic indices than the typical 
SFGs (see also Figure [S]) suggesting that they have ex¬ 
perienced a substantial structural transformation. To 
account for the presence of these compact population on 
the Ee and Ei relations, we fit the logE — logM* distri¬ 
bution using a second order polynomial, which allows for 
a change in the slope at the high mass end. The thin blue 
lines in Figure [5] indicate that the relations become in¬ 
deed steeper at log(M/M 0 )> 10.5 in both Eg and Ei (see 
also Table [U due to the increasing number of SFGs with 
higher surface densities and Sersic index (green circles). 

In the following, we refer to any evolutionary paths 
upwards from the E-MS due to phase(s) of high-efficiency 
core growth as compaction track(s). We emphasize that 
compaction is, primarily, a core building process, and 
thus it is more efficient at increasing Ei than Ee, as the 
latter depends also on the overall galaxy size. See for 
example the SFGs with dense cores selected within the 
scatter of E)^ (green circles in Figure H]) which have a 
larger scatter in Eg. We further discuss this difference in 
§3.5. 

3.4.2. Trends in hydrodynamical simulations 
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In order to provide further support for the structural 
evolutionary paths discussed above, and obtain some in¬ 
sight on the physical mechanism(s) that cause them, 
we study the evolutionary paths of a sample of high- 
resolution hydrodynamical simulations in the same pa¬ 
rameter space. Figure [5] shows the evolutionary tracks 
in Si and Se vs. mass for th e Vela set of 35 simulations 
descri bed in iCeverino et al.l (|2014ll and iZolotov et al.l 
(|2015ll . In the latter, the authors analyzed the SFR, 
structure and kinematics of the simulated SFGs, show¬ 
ing that the driving force behind the most significant 
changes in all these properties is a phase (or phases) 
of dissipational “compaction” caused by, e.g., mergers, 
disk instabilities, interactions, counter rotating streams 
or tidal-compressions, which trigger intense episodes of 
gas inflow and SFR, as ment ioned in the pr e vious section. 

Following the approach in IZolotov et al.l (I2015f l. we di¬ 
vide the sample in two groups according to their mass (in 
Figure 6, high-mass in black, low-mass in grey), and we 
show their median evolutionary tracks at early (dashed; 
z > 2.5) and late (solid; z < 2.5) times. The late evo¬ 
lution matches the redshift range of the galaxy sample, 
shown as blue and red contours. Qualitatively, the tracks 
have a similar behavior in both panels, and match the ob¬ 
served distributions. To first order, the overall evolution¬ 
ary paths can be described in terms of two phases: an 
early phase of steady structural growth approximately 
following the slope of the E-MS (a ^ 0.6, /3 ~ 0.9), 
and a compaction phase, which causes a steep increase 
in Seand Ei (a ~ 1.5, /3 1.3). The steep increase 

in Ee also implies a size shrinkage, parti cularly in the 
high- mass simulations (see also Figure 9 in IZolotov et al.l 
1201511 . Overall, the dependence on re, which can change 
rapidly due to gas accretion and/or interactions, implies 
a larger spread in Ee with respect to Ei for a given stellar 
mass. 

The main differences in the black and gray tracks are: 
1) a downsizing effect, i.e., the most massive galaxies 
evolve earlier and faster; 2) the evolutionary tracks of 
massive galaxies have a higher normalization; 3) the 
main compaction event is stronger and happens earlier 
for the most massive galaxies. The decline in the nor¬ 
malization is the result of a decline in the gas fractions , 
and thus SFRs, as a function of time (iDekel et al.ll2013l : 
IZolotov et al.l[201^ . 

The different loci of SFGs and quiescent galaxies (blue 
and red contours) also results from the gradual decline 
in the gas reservoirs and the differential evolution with 
stellar mass. The quiescent locus consist, almost ex¬ 
clusively, of massive galaxies that experienced a strong 
compaction event at z > 2.5 (black lines in Figure |6]). 
These galaxies quenched star-formation in the core, but 
continued to grow i n stellar mass due t o star-formation 
in a re-grown disk (IZolotov et al.l[2(TT5ll . The new disk 
causes a sudden increase in size, which leads to decreas¬ 
ing tracks with mass in Ee, but nearly flat tracks in Ei. 
The star-forming locus consist mostly of low mass galax¬ 
ies which populate the high-mass region at lower redshifts 
z < 2.5 (gray lines). These galaxies experience weaker 
compaction events, after which the core resumes star- 
formation and continues on an evolutionary path with a 
similar slope as the Ei-MS. Later, these SFGs can have a 
secondary compaction event and arrive on the quiescent 
sequence from below, as illustrated by the green track of 


Velal3. This galaxy has the latest compaction event and 
thus spends most of its life in the E-MS. Interestingly, 
during the E-MS phase, the gas mass in the central 1 kpc 
remains almost constant, as expected in the case of a sim - 
ple “bathtub” model evolution fe.g.. IDekel et al.ll2013jl . 
This “bathtub” phase strengthens the notion that the 
structural main sequence describes a phase of smooth, 
steady-state evolution that coincides with the SFR main 
sequence. 

Lastly, note that despite the overall agreement in ob¬ 
servations and simulations, there are some differences 
between the simulated tracks and the observed param¬ 
eter space in E vs. mass, particularly at high-z. At 
z > 2.5, all massive SFGs evolve directly into the high- 
density, quiescent region (red contour) due to strong 
compaction events, which also cause significant mass 
growth (AlogM* > 1 dex). This evolution under¬ 
predicts the number of massive SFGs (log(M/MQ)> 10) 
with low Ee and Ei at z > 2.5, which suggests that the 
E-MS evolutionary channel is underrepresented, or hap¬ 
pens only at low redshifts, in this sample simulations. 


3.4.3. Summary of evolutionary paths 

Based on the excellent agreement in the observed 
galaxy distributions and model evolutionary tracks in the 
log E — log M* diagrams, we speculate that the overall 
long-term structural evolution of SFGs can be expressed 
in terms of 2 fundamental phases, namely, a E-MS and 
compaction. These phases approximately follow linear 
tracks in logEi^e cx; [a,/3] log M*. The E main sequence 
{a 0.6, l3 0.9), is a relatively smooth phase of struc¬ 
tural growth, consistent with a period of steady-state gas 
accretion and galactic scale star-formation. As a result, 
the central and effective densities increase with time, as 
does the galaxy size, i.e., similar to the typical inside-out 
growth of a disk. Compaction, is a phase (or phases) of 
enhanced core-growth {a > 1, /3 > 1) fueled by strong 
gas infall to the galaxy center, as a result of gravita¬ 
tional instabilities. In this phase, the central and effec¬ 
tive densities increase steeply, and the galaxy’s effective 
radius can decrease due to the large increase in stellar 
mass close to the center or due to a structural collapse. 
Overall, the intensity and duration of compaction events 
decreases with time due to the similarly decreasing gas 
fractions. 

For quiescent galaxies, the simulations reproduce 
the expected evolutionary path of decreasing Ee and 
relatively constant Ei with increasing mass. How¬ 
ever, this path is fueled by low levels of in-situ star- 
form ation, which appears to co nt radict the fast quench¬ 
ing (Ivan de San de et al.1 120131 : iBezanson et al.l 120131 : 
lOnodera et al.ll20 14 ) and the lack o f diffus e star-forming 
components (iSzomoru et ^ I2011L l2012fl in quiescent 
galaxies at z > 1.5. The leading theory is that those 
quiescent galaxies grow due to minor mergers in absence 
of further star-formation. Alternatively, it could be that 
some of the SFGs with high surface densities are reju¬ 
venated quiescent galaxies with re-grown star-forming 
disks. Nonetheless, given the rapidly increasing number 
of quiescent galaxies (see § 13.5.311 and the lack of appar¬ 
ent disk signatures in those SFGs (see i) l3.6p . rejuvenated 
quiescent galaxies with extended star-forming disks can 
only make a small fraction of the sample. 
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Fig. 7. — ASFR vs. AE plots from the SFR-MS and the Se (top) and Si (bottom) quiescent structural relations as a function of redshift 
for galaxies with log(M/M 0 ) > 10.3. The x-axis is normalized to the zero-point of the quiescent relations to illustrate the evolution with 
time. Blue and green circles indicate SFGs outside and within the 2cr scatter of S^ (dashed lines). We refer to the latter as compact SFGs. 
The relative fractions of SFGs and quiescent galaxies found within the 2a scatter of E^ and S^ (dashed lines) are indicated in the bottom 
left. The L-shaped nature of the distributions indicate that SFGs become compact before they quench their star-formation. Moreover, 
the L-shape persist across redshifts, implying a universal process of compaction followed by quenching (see Figure ( 8 ) 1 . The distributions 
in the top and bottom panels are very similar; however, EJ* exhibits a tighter scatter and a much slower evolution of the zero-point. The 
purple line shows the compactness threshold in E^from B13 which is only efficient at z > 2. A fraction of the compact SFGs selected with 
AEJ are not compact in EJ. This indicates that those SFGs have compact cores, but have larger r^ than typical quiescent galaxies. Thus, 
a threshold in E J is a more efficient selection criterion to identify compact SFGs and quiescent galaxies. 


3.5. Relative distances from the SFR-MS and the 
quiescent structural relations: compaction and 
quenching 

In this section we study the relative distributions of 
galaxies with respect to the SFR-MS and the struc¬ 
tural relations to identify candidate quenching galaxies 
as a function of redshift. Qualitatively, this analysis is 
similar to studying the distribution with respect to the 
SFR- and E- main sequences. However, using the E^^ 
frames the selection around quiescent galaxies, searching 
for SFGs with similar structural properties to them. 

3.5.1. Compact SFGs as progenitors of quiescent galaxies 

As discussed in the previous section, the higher nor¬ 
malization of the quiescent E^ and E^ relations with 
respect to those of SFGs indicates that quenching is pre¬ 
ceded by an increase in the surface density above a cer¬ 
tain threshold. However, as shown in Figure[2l such char¬ 
acteristic density scales with stellar mass, challenging the 
simple notion of a fixed quenching thr eshold at a given 
surface density or stel lar mass (e.g., iKauffmann et al.l 
120031: iFranx et al.ll2008f ). In turn, the most effective qui¬ 
escent criterion is a selection with respect to the E^ or 
E^ sequences with stellar mass, i.e., a relative offset from 
the structural relations. Such relative selection includes 
fewer SFGs. However, there is always overlap with quies¬ 
cent galaxies at every redshift. A possible interpretation, 
outlined in the previous section, is that those overlapping 


SFGs acquire quiescent morphologies while they are still 
star-forming as a result of a compaction process, i.e., a 
structural transformation that increases the central den¬ 
sity and Sersic index, and reduces the size before quench¬ 
ing star-for mation. Such an evolution ary sequence was 
confirmed in iBarro et al.l (|2013l . [^14al hereafter B13) for 
SFGs at z ~ 2. In BI3, the authors used a selection on 
specific SFR and relative distance to the quiescent size- 
mass relation to identify compact SFGs at z ~ 2, finding 
that those galaxies have similar sizes, Sersic indices and 
spheroidal morphologies as the quiescent population. 

Building on this idea, in this work we define com¬ 
pact SFGs as those SFGs (ASFRms > —0.7 dex) found 
within the ^ 2x the scatter of the quiescent struc¬ 
tural scaling relations at a given redshift, AE^j^ = 
log Eg , 1 — logE^;^(z) > —2CT(logE‘3), where Eg^i is ei¬ 
ther the central or effective mass density, and we use 
2cr(logE^) = 0.5 dex and 2cr(logE^) = 0.2 dex, re¬ 
spectively. This definition differs from previous works 
where compact is an absolute term to iden tify the small¬ 
est galaxies at high-z (cg < I kpc; e.g., iGassata et al.l 
1201 IL [20T^ . Here, compact is a relative term referring 
to the densest/smallest SFGs at every redshift. Panels A 
and B of Figure[7]illustrate the selection of compact SFGs 
in Eg and Ei (dashed line). In the x-axis we add the 
zero-point at each redshift to illustrate the different time 
evolution in the normalization of these relations. Panel 
A shows also the compactness threshold of BI3 (purple 
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Compaction 



7 8 „ 9 JO 11 

AE°+10[M^kpc-2] 


Quenching Sequence 



8.5 9 „ 9.5 ,10 10.5 

AE°+ 10[M^kpc“1 


1SFR & £ Main Seq. 

ASe > 0 ; ASi > 0 
Are > 0; An ~ 0 

2. - Compaction 

ASe > 0; ASi > 0 
Are < 0; An :§> 0 

3. - Quenching 

ASe > 0; ASi ~ 0 
Are < 0; An ~ 0 

4. - Passive core 

ASe < 0; ASi ~ 0 
Are ^ 0; An ~ 0 


Fig. 8. — Stacked ASFR-MS vs. ASe (left) and ASj (right) for all massive galaxies at 0.5 < 2 : < 3, color coded by Sersic index. The 
AS values are normalized at Sj=S^=10. The gray lines show the best-fit linear relations for SFGs and quiescent galaxies. The L-shape 
and the intrinsic scatter of these diagrams is independent of redshift, which suggests that all massive galaxies follow the same universal 
evolutionary path from star-forming to quiescent. The numbers and the legend on the right side summarize the main phases of such 
evolution: 1) The SFR and S- MS phase is a nearly horizontal track in which galaxies increase both their core mass and effective radius 
at relatively constant Sersic (i.e., a disk growth); 2) Compaction is a phase of enhanced core growth, that increases the Sersic index. Both 
1) and 2) follow horizontal tracks, however, only compaction causes a substantial increase in n, pushing SFGs galaxies towards the knee of 
the relation. 3) Quenching is a nearly vertical track that indicates that the full shut down of star-formation happens at maximum Ee and 
El. 4) In the passive core phase Ei remains constant (or declines slowly due to dying stars), while Eg decreases faster due to size growth 
or new arrivals, causing a large evolution in the normalization of EJ (Figure 0. The slight tilt of the E-MS to compaction phase indicates 
that the onset of quenching starts begins before reaching the maximum stellar density. 


line), which, by definition, matches < 0.5 dex at 

z ^ 1.8. However, as the normalization of S^declines 
with time, a fixed selection threshold gets progressively 
fewer galaxies. In turn, ASFRms — AE^ is essentially a 
redshift-independent extension of the method to identify 
star-forming progenitors of quiescent galaxies. 

The relative distributions of SFGs and quiescent galax¬ 
ies in both panels of Figure [7] are very similar. By def¬ 
inition, both AE^ and AE^ select all quiescent galax¬ 
ies within the typical scatter of the structural relations. 
A Ee selection exhibits a few more catastrophic out¬ 
liers, most likely due to extr eme values of the Sersi c 
also affecting the Cg (see e.g., Ivan der Wei et al]l2012h . 
Nonetheless, both Ee and Ei identify the bulk of the 
quiescent population. The selection of compact SFGs is 
also largely consistent. However, Ei selects ^ 10 — 15% 
more compact galaxies at each redshift. These compact 
SFGs in El scattered outside the Ee selection have both 
dense centers and extended star-forming profiles, which 
leads an overall higher in a single Sersic fit (e.g., 
iBruce et al.ll2012l) . Non etheless, as discussed for exam- 
ple in iFang et all (|201,'1[1 . these galaxies will eventually 
increase their Ee as the extended star-forming disk fades. 
Thus, they are bonafide quiescent progenitors. The lower 
efficiency of Ee for selecting compact SFGs is a result of 
Ce not being monotonic with time (i.e., Vg both increases 
and decreases with redshift). 

In summary, Ei is closer to a clock (it only increases), 
and it also exhibits a tighter scatter and a slower evo¬ 
lution the normalization, so we conclude that AE)^ is 
a more efficient criterion to identify compact SFGs. In 
fact, as shown in panel B of Figure 0 a single threshold 
of El — 0.65 log(M* — 10.5) > 9.5 M 0 /kpc^ identifies the 


majority of compact SFGs and quiescent galaxies at red- 
shifts z > 0.5. Hereafter, we will refer to compact SFGs 
as those selected in Ei. 

3.5.2. A universal compaction-quenching sequence 

The remarkable similarity in the galaxy distributions 
of Figure [7] as a function of time suggest that the relative 
distance from the SFR-MS and E^,^ relations describes 
a universal evolutionary sequence for massive galaxies, 
which is independent of redshift. Figure [ 8 ] illustrates this 
sequence showing the stacked distribution in ASFRms 
vs. AE^ for the redshift range 0.5 < z < 3.0. The 
upside-down, L-shaped sequence summarizes the notion 
that forming a dense ste llar core is a pre-requisite for 
quen c hing star-formati on (jGheung et al.ll20i^ I Bell et aP 
120121 : iFang et HI 120131 1 . The color code in Sersic index 
emphasizes that the evolution requires both the growth of 
a dense stellar core, and a structural transformation from 
an exponential (disk) profile to more a concentrated (n > 
2) Sersic profile. Incidentally, Figure | 8 ] also illustrates 
why the Sersic index is a better quiescent indicator than 
a con stant threshold in Eg, Ei or velocity dispersion, cr 
(e.g., iBell et aDl2012[ l. 

In terms of the galaxy evolutionary paths discussed 
in §E11 the two branches of the L-shaped sequence de¬ 
scribe two fundamental transitions from the E- and SFR- 
main sequences, respectively: compaction and quench¬ 
ing. The horizontal branch of SFGs represents constant 
evolution along the SFR-MS, but an eventual departure 
from the E-MS due to a compaction event. The E-MS 
describes the growth of an exponential disk (n ^ 1 ) due 
to in-situ star-formation, which increases Ei, Eg and 
the overall size of the galaxy. Compaction involves a 
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stronger increase in Si and Se, leading to a steeper mass 
profile (n > 2) and a smaller r^. The vertical branch 
shows that quenching (defined as ASFRms > —0.7 dex) 
takes places at maximum central density in compact 
SFGs with similar morphologies to the quiescent pop¬ 
ulation, i.e., compaction precedes the full shut down of 
star-formation. The slight tilt in the horizontal branch, 

ASFRms-0.20AS?, ASFRms -0.25ASQ (gray 

lines), however, suggests that while quenching is com¬ 
pleted at maximum central density it starts some time 
before. 

Following the interpretation of B13, we characterize 
the evolutionary pace along compaction-quenching path 
as declining from fast to slow as a function of time. The 
gradient in quenching speed a rises naturally from the de¬ 
cline in sSFR with time (e.g.. iSoeagle et al]l2014ll . Since 
A log Ei^e oc A log M oc sSFR(z), the central densities in¬ 
crease more slowly with time. The la rger SFR at high-z is 
associated with high er gas fractions (jTacconi et al.ll2n'T(Tl : 
iTacconi et al.l[2M^ . which are likely related with larger 
accre t ion rates from the dark matter haloes (|Dekel et al.l 
120131 : iDekel fc BurkertI 1201411 . Simple stability argu¬ 
ments predict that such gas rich galaxies are prone also to 
stronger gravitational instabilities, and thus favor more 
expedite compaction and quenching processes. In such 
case, the tilt in Figure HI could be due to gas starva¬ 
tion after an instability-induced starburst, i.e., a wet- 
inflow that causes a peak in SFR at maximum gas den¬ 
sity, and declines progressively with the gas supply, while 
the core mass incre ases le.g.. iFeldmann fc Maveil [20151 : 
iZolotov et al.ll201^ . The truncation of the gas supply 
to the galaxy center can be caused by a combination of 
feedback (e.g., outflows) and stab ilization of the disk due 
to the increasing central density (iMartig et al.l 1200^ . 

At lower redshift, the smaller gas fractions and 
longer dynamical timescales suggest that compaction 
mechanisms become a mixture of weaker instabilities, 
which can still cause enhanced ga s inflows or inward 
migration of stellar cl umps (e.g., Dekel et ^ l2009bl : 
iBournaud et al.l 120111: iGenel et alj_ 201211 . and secular 
processes (iKormendv fc Kennicuttll2004[ l associated with 
torques and dynamical friction in the presence of bars 
and spiral arms. The latter play an important role in¬ 
creasing the central density in SFGs that already ex¬ 
hibit relatively quiescent centers, i.e., lacking enough 
star-format ion to sustain the core growth requi r ed to 
reach S? ( Wuvt s et al.l I2012L l2013t iLang et Si 12014 
iBruce et al.il2014bi ah. In those SFGs, quenching is also 
expected to be a slow process, related with ga s consump¬ 
tion in the star-forming d isk (fading; e.g., iFang et H] 
120131 : iTacchella et al.ll^OlSH . The slow quenching process 
depends also on additional mechanisms to prevent fur¬ 
ther gas accretion into the disk, e.g., virial shock heating 
in ma ssive haloes (|Groton et al.l[2r)0ft IDekel fc BirnboimI 
1200011 or AGN feedback (. As mentioned in the previous 
section, some of these SFGs are compact in Ei but not 
in Ee. Therefore, the tilt in ASFRms vs. AE^ can be 
caused by the size shrinkage (rg) associated with fading 
rather than with an increase in their core mass. 

3.5.3. Number density of compact SFGs and quiescent 

galaxies 

The existence of a universal compaction-quenching di¬ 
agram and the increasing number density of quiescent 



Fig. 9.— Redshift evolution of the number density of massive 
(log(M/MQ)> 10.3) compact SFGs and quiescent galaxies. Com¬ 
pact SFGs are selected using the criterion. The red triangles 
indicate the evolution in the nu mber density o f quie scent galax¬ 
ies in different mass ranges from lMuzzin et al.l II2013I) . The green 
and red lines show the best-fit evolutionary model to the observed 
number of compact SFGs and quiescent galaxies. The model is 
based on the assumption that quiescent galaxies are descendants 
of compact SFGs that have a characteristic quenching timescale 
aQ, but can also increase in number due to compaction events in 
more extended SFGs that have a frequency This simple 

model that assumes a continuous replenishment of compact SFGs 
at every redshift can account for the observed evolution in the two 
populations. 


galaxies as a function of time supports the idea that the 
L-shaped diagram is indeed an evolutionary sequence, 
i.e., compact SFGs in the knee of the relation are im¬ 
mediate progenitors of the quiescent galaxies at lower 
redshifts. Figure H] shows the evolution in the number 
density of compact SFGs, selected using AE^ > —0.2, 
and all quiescent galaxies, regardless of their structural 
properties, since z ~ 3. The latter grows approximately 
as a power-law, ^ io-0 -5(i-bz)^ good agr eement 
with previous results (e.g., iMuzzin et al.l 1201311 . As¬ 
suming that quiescent galaxies are descendants of com¬ 
pact SFGs at higher redshift, the number of quiescent 
galaxies is the cumulative distribution of quenching com¬ 
pact SFGs as a function of time. To quantify the 
relative numbers of these two populations we assume 
that the number density of compact SFGs also follows 
a power-law evolution as the result of 2 opposite pro¬ 
cesses that cause either an increase (due to compaction) 
or decrease (due to quenching) in their total number: 
dNsFc/dz = AcA^sfg + AqAsfg- In this model, the 
evolution in the number of quiescent galaxies is propor¬ 
tional to the number of quenching SFGs, and thus in¬ 
herit the power-law dependence: dN^/dz = — A^A^seg 
—T here are two additional parame¬ 
ters which determine the initial number compact SFGs 
and quiescent galaxies, Nfp® and N^q. We set nJq = 0 at 
zO = 4, and we fit the other 3 parameters to the observed 
number density of compact SFGs and quiescent galaxies, 
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Fig. 10. — Core to total mass ratio (C/T) vs. stellar mass for a representative sub-sample of massive galaxies at different redshifts. The 
color composite ACS/WFC3 zJH images illustrate the evolution in the morphology of SFGs and quiescent galaxies as a function of time. 
The color of the border around the postage stamps indicate the location of the galaxy in the three regions of the compaction-quenching 
diagram (Figure [S}. The solid blue and red lines indicate the best fit Ei relations for SFG and quiescent galaxies at each redshift. The 
dashed lines indicate the typical width of the relation, which is the threshold used to select compact SFGs. The insets show the same 
diagram color coded by Sersic index for all massive SFGs. SFGs evolving in the E-MS decrease their C/T with time as the disk component 
becomes more prominent (i.e., the disk grows faster than the core). Similarly, quiescent galaxies decrease their C/T with time as they 
grow larger stellar envelopes. However, SFGs that have a compaction event, experience a significant core growth relative to the total mass 
(increasing their C/T), moving upwards from the E-MS to the quiescent relation. The visual appearances of compact SFGs (g reen ) are 
very similar to those of quiescent galaxies at all redshifts, but are most different from other SFGs only at high-z (see also Figure [TTI l. 


which yields \i = 1.27 ± 0.15, \d = —0.75 ± 0.12 and 
= 5-10“® Mpc“^. Thus, the characteristic quench¬ 
ing timescale for a compact SFG in units of redshift is 
^ 1 / 2 ,Q = In(lO)/ln(2)/AQ = 0.38, which implies that the 
quenching time increases from tq = 700 Myr to 1.1 Gyr 
at z = 3 and z = 1. In spite of its simplicity, this model 
provides a more realistic approximation than previous 
assumptions that all compact SFGs form and quench in 
discrete intervals of time. In this case, the number densi¬ 
ties vary continuously while preserving the evolutionary 
connection as a characteristic quenching time. 

3.6. Galaxy morphologies in the universal 
compaction-quenching sequence 

As discussed in § 13.5[ the shape of the compaction¬ 
quenching sequence is independent of redshift. However, 
the evolutionary pace along the sequence declines with 
time as compaction, and most likely quenching, processes 


become slower. Furthermore, the evolution in the nor¬ 
malization of Se clearly indicates that the structural and 
morphological properties of galaxies in the sequence de¬ 
pend on the redshift. 

Figure [TU] illustrates the evolution in the visual appear¬ 
ances of SFGs and quiescent galaxies as a function of 
time in the core-to-total vs. mass diagram. The thumb¬ 
nail frames are color coded according to the location of 
the galaxies in the 3 regions of the compaction-quenching 
diagram (Figure El). Overall, the evolutionary tracks in 
this diagram follow decreasing trends in G/T for both 
SFGs (i.e., inside-out growth of a disk), and quiescent 
galaxies (i.e., stellar halo growth). Only SFGs experi¬ 
encing a compaction event (from blue to green) follow a 
trend of increasing G/T. Note that, conceptually, C/T 
is similar to a bulge-to-tota\ ratio. However, bulges can 
typically grow beyond a radius of 1 kpc, in fact, recent 
results indicate that both the bulge size and its mass 
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Fig. 11. — Postage stamps in the ASFR—AS^ plane, showing galaxies in three different phases of the compaction-quenching sequence 
(see Figure 0: extended SFGs (blue), compact SFGs (green) and quiescent galaxies (red). Each panel shows galaxies with increasing 
stellar masses as a function of time (left to right) following the expected mass growth from z = 3 to z = 0.5 predicted by semi-analytic 
models . The numbers on the inset indicate the median n and re of the galaxies in each group. The evolutionary sequence, discussed in 
^ 13.5.21 suggests that some SFGs grow along the S-MS from z = 3 to z = 0.5 (blue), while others experience a compaction event (green), 
departing from the S-MS. Gompaction produces a steeper mass profile, resulting in a strong increase in Ei and Eg. After compaction, 
galaxies quench (red), showing only minor changes in their structural and visual appearances relative to their compact SFG progenitors. 
The E-MS describes the growth of a disk, starting from a small, irregular morphology and becoming progressively more settled with time, 
increasing both its size and central density (bulge growth). The pre- and post- compaction appearances of SFGs at high-z are radically 
different, suggesting the action of a strongly dissipative process. However, at low-z, many compact SFGs exhibit disks, indicating that 
compaction affects mostly the central region of the galaxy. The overall appearance of quiescent galaxies evolves from “naked” to “clothed” 
cores (i.e., larger Ve at similar Ei). Such evolution is likely driven by size growth in quiescent galaxies formed at high-z, and the formation 
of new quiescent galaxies from compact SFGs with larger sizes. 


correlate with the total mass of the galaxy, i .e., B/T in¬ 
creas es with mass, as opposed to C/T (e.g., iLang et HI 
I2014f) . Qualitatively, the visual appearances agree with 
this evolution. The bulk of SFGs exhibit larger are more- 
settled disks as a function of time, while also increasing 
their central densities. The steeper slope of E®'’ seems 
clear in this diagram (dashed blue line), particularly at 
z < I, where slow compaction processes push most SFGs 
above C/T> 10%. Quiescent galaxies also show larger 
{puffed up) appearances with time with values declining 
from C/T ~ 60% at z = 2.6 (“naked”-cores) to 35% and 
z = 0.75 (“clothed”-cores). 

To further illustrate the change in visual appearances 
within an evolutionary path. Figure [TT] shows images of 
galaxies with increasing stellar mass a function of time in 
the 3 regions of the compaction-quenching sequence. The 
increase in stellar mass with redshift follows the average 
mass growth in galaxies with log(M/MQ) = 10.5 at z = 
2.5. There several methods to estim ate the typical stel¬ 
lar mass growth with redshift (e.g. Ivan Dokkum et al.l 
120101 : iPatel et al.ll2013l : [Marches iii et al.||201^. Her e we 
adopt the evolutionary tracks of iMoster et al.l (j2013l l de¬ 
termined from semi-analytic models, with average mass 
growth of AM* = 0.7 dex from z = 2.5 to z = 0.75. Since 
galaxies of a given mass can be in any of the three stages 
at a given redshift, we assume that at each redshift some 
SFGs remain in the left region (blue) growing along the 
Ei-MS, others move rightward due to compaction (blue 
to green), and some compact SFGs quench downward 
after reaching a maximum central density (green to red). 

The morphologies and visual appearances change 
strongly as a function of redshift: 1) at z > 1.4, com¬ 
paction involves a significant structural transformation 
going from irregular, clumpy SFGs to compact “naked” 
spheroids. Compact SFGs quench having almost iden¬ 
tical morpholo gies than quiescent galaxies at those red- 
shifts (see also iBarro et al.ir2014af) . Extended SFGs in 
the E-MS (blue) increase their size, central density and 
total mass with time. However, as total mass grows 
faster than the core mass, their C/T decreases. By z < 2 


these galaxies exhibit clear disk-like structures and dis¬ 
tinct cores with redder colors (see alsoiWuvts et al.ir2ni2l : 
ILang- et H] 12014 iTacchella et al.ll2ni5h . 21 at z < 1.4, 
extended SFGs show larger, more settled disks and sim¬ 
ilarly more massive cores. Compact SFGs and quiescent 
galaxies also have denser cores, but, in contrast with 
the high-z appearances, they exhibit diffuse extended 
components (some clearly star-forming). This suggests 
that the compaction process at lower redshift preserves 
the already existing disk structure. The evolution to¬ 
wards disk-like morphologies with progressively bigger 
cores eventually results in the prominence of disks with 
bulges among massive SFGs and qu iescent galaxies at 
,z 1 reported in p r evious works (e.g..lBundv et al.ll2ni(1l: 
Bruce et al.l l20i a iBruce et al.1 I2n i4a: iBuitrago et al.l 
201,31 : iMcLure et al.1 1201,31 : ILang et al.l 120141: Huertas- 
Company et al 2015). As noted in previous works, how¬ 
ever, the progenitor bias makes it nearly impossible to 
distinguish between puffed-up, older galaxies and newly 
quenched galaxies with larger sizes attending only to 
their morphologies. 

4. SUMMARY 

We analyze the star-formation and structural prop¬ 
erties of massive galaxies in the CANDELS/GOODS-S 
field to study the relation between stellar mass and struc¬ 
tural growth, and the role of the latter in the quenching 
of star-formation since z ~ 3. We characterize the struc¬ 
tural properties as a function of redshift by studying the 
correlations in the mass surface density within the ef¬ 
fective radius, Ee, and within the central 1 kpc, Ei, vs. 
stellar mass. Ei traces the stellar mass growth in the 
galaxy core, and thus it is close to the concept of a cos¬ 
mic clock (i.e., it increases with time). Ee , however, de¬ 
pends on the relative balance between stellar mass and 
size growth, and thus it exhibits positive and negative 
fluctuations. 

We find that SFGs and quiescent galaxies follow clear 
and distinct correlations in Ee and Eivs. stellar mass 
since z ~ 3. These correlations are well-described by lin- 
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ear relations in log-log space. The slopes and scatter of 
these relations are relatively constant with time, while 
their normalizations decline (see Table [T|) . The scatter 
in the Ei structural relations is ~2x tighter than in the 
Ee relations for for both SFGs and quiescent galaxies. 
For SFGs, the normalizations in Ee and Ei decrease by 
less than factor of ^2 from z = 3 to z = 0.5. For qui¬ 
escent galaxies, the decline in Ee is 5x larger than in 
El (~ 0.3 dex vs. ~ 1 dex, respectively). The differen¬ 
tial evolution of the normalization in E^ and E^ is in¬ 
consistent with a simple minor merger scenario in which 
the core mass remains constant while the size increases. 
However, the normalizations at redshifts z = 3to z = 0.5 
agree well with one another if the mass profiles of quies¬ 
cent galaxies follow a single Sersic with n ~ 4. Thus the 
large increase with time in the size of quiescent galaxies 
is consistent with just modest decline in Ei at constant 
Sersic. Such decline can be caused by the formation of 
new quiescent galaxies with lower density cores or mass 
loss due to passive evolution in the already existing pop¬ 
ulation. 

Based on the slow decline in the normalizations of the 
structural relations for SFGs, we speculate that these 
galaxies follow evolutionary paths along the Egq corre¬ 
lations. We define these paths as the E main sequence, 
following the paradigm of the SFR-MS as a smooth phase 
of stellar and structura l growth (e.g., lElbaz et '311(20111 : 
iRodighiero et al.]l2011f) . The evolution in the E-MS is 
consistent with the inside-out growth of an exponential 
disk due to in-situ star-formation, i.e., an increase in both 
the core mass and the overall size of the galaxy. 

At every mass and redshift, quiescent galaxies have 
steeper mass profiles (higher Sersic) and higher surface 
densities than SFGs. This implies that growing a dense 
stellar co re is a pre-requisite for quenching star-formation 
(see a l so.lGheung et al. 11201^ I Bell et al.ll201^ iFang et HI 
[Ml I van Dokkum et al.l 120141) ' Thus, the immediate 
star-forming progenitors of quiescent galaxies must ex¬ 
perience a phase of stronger core-growth, relative to 
the E-MS. We define this phase(s) of fast increase in 
El, Ee and n as compaction. The compaction phase is 
typically associated with dissipational processes ranging 
from major mergers (iHopkins et al.l[2006t I Hopkins et al. 
200^ to violent gravita t ional in stabilities (iDekel et al. 
2009ai: IDekel fc BurkertI 12014 iCeverino et al.l 1201011 . 
and weak er secular instabilities ( e.g., bars and spi¬ 
ral arms; iKormendv fc KennicuttI 120041) . The evolu¬ 
tionary tracks of massive SFGs in recent hydrodynam- 
ical simulations exhibit an exce llent agreement with 
the E-MS/compac tion scenario (ICeverino et al.l 120151 : 
iZolotov et al.ll20ll) . The simulations suggest that SFGs 
depart from a steady-state evolution along the E-MS as 
a result of dissipational compaction events triggered by 
intense episodes of gas inflow. The strength and duration 
of these events declines wit h time following a decline in 
the gas reservoirs of SFGs (|Dekel et 31.1120131) . 

We find that the 2D distribution of SFGs and quies¬ 
cent galaxies relative to the SFR-MS and the quiescent 
structural relations, ASFRms — exhibits a uni¬ 

versal L-shape that is independent of redshift (Figure [7]). 
Each branch of this distribution describes a fundamen¬ 
tal transition in the evolution of SFGs, namely: com¬ 
paction and quenching. The horizontal branch describes 
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Fig. 12. — Structural evolutionary path of massive galaxies since 
2 3 in terms of their core (r < 1 kpc) to total mass (C/T) vs. 

stellar mass. The typical SFG follows a structural “main sequence” 
that describes the inside-out growth of a disk (blue-framed images), 
i.e., increasing its stellar mass and Ve as does its central density 
(bulge growth). Prior to quenching, SFGs experience a dissipative 
compaction phase of enhanced core growth (green-framed images) 
fueled by the inward migration of gas into the central region. Com¬ 
paction events are strongest at high-z, due to higher gas fractions, 
forming extremely compact SFGs on very short timescales. At 
low-z the strength of compaction diminishes, the core growth be¬ 
comes slower and it preserves the underlying disk structure. Newly 
quenched galaxies at each redshift (red-framed images) are direct 
descendants of compact SFGs, which shut down star-formation af¬ 
ter forming a dense stellar core with no further structural changes. 
Quiescent galaxies grow an extended stellar halo due to minor 
merging, decreasing their C/T. Quenching timescales depend on 
the strenght of the compaction process, which drives a major gas 
consumption event. Thus, qualitatively, quenching follows the fast 
shut-down of a compact core at high-z and the slower fading of a 
star-forming disk at low-z. 


the growth in size and core mass of a SFG along the E- 
MS, followed by a compaction event that causes a strong 
increase in surface density and Sersic index. Gompaction 
marks the departure from the E-MS towards the knee of 
the L-shaped distribution. The vertical branch describes 
the quenching of star-formation at maximum core and 
effective surface density. The latter implies that the for¬ 
mation of a dense core precedes the quenching of star- 
formation. Therefore, compact SFGs in the knee of the L- 
diagram at every redshift are the most likely progenitors 
quiescent galaxies at lower-z. Owing to the mild decline 
in the normalization of the E^ relation with time, a con¬ 
stant threshold of E^ = Ei — 0.651og(M* — 10.5) > 9.5 
identifies compact SFGs at any redshift z > 0.5. 

Although the L-shape of the compaction-quenching se¬ 
quence is independent of redshift, the rate of growth 
along the sequence and the strenght of the transforma¬ 
tion processes declines with time following the similar 
decline in the over all SFR (e.g., iWhitaker et al.l l20l4 
iSneagle et al.ll2?)T^ . As a result, the morphologies of 
galaxies in each region of the sequence are substantially 
different with time. Figure [T^] describes qualitatively 
the change in morphologies and core-to-total ratios for 
galaxies following fast-to-slow evolutionary paths in the 
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compaction-quenching sequence as a function of time. 

At high-z, SFGs have larger gas reservoirs and 
thus follow a faster mass growth (iTacconi et al.l 120101 : 
iTacconi et al.l 1201, *111 . Those SFGs are also prone 
to stronger gravitational instabilities and expedited 
compaction processes capable of transforming SFGs 
with clumpy, irregular morphologies into compact 
SFGs with spheroid-like morph ologies, high Sers i c and 
and high core-to-total ratios (iDekel fc BurkertI [2014 
Wellons et al.ll201^ iGeverino et al.ll201,^ : I Zolotov et al l 
201,^ . l.e., “naked” star-forming cores, which appear 
to quench fast turning into compact quiesc e nt galaxies 
(iBezanson et al.l 120131 lOnodera et al.l 12014 iBelli et ^ 
l2014bHBarro et al.ll2015ir with similar masses and struc¬ 
tural properties. 

At low-z, SFGs have lower gas fractions and SFRs, 
which implies a slower mass growth and weaker gravi¬ 
tational instabilities. This is consistent with SFGs hav¬ 
ing more settled, disk-like morphologi es, larger sizes, and 
progressively denser cores (see also, jBruce et al.l 120121 : 
iLang et all 12014 iTacchella et al.ll2015h . As opposed to 
high-z, the similar morphologies of compact SFGs (green¬ 
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framed images) and SFGs in the E-MS (blue-framed im¬ 
ages) suggests that compaction at low-z causes an en¬ 
hanced core growth, increasing G/T within an already 
existing disk. Many SFGs at low-z exhibit blue disks 
and red cores, suggesting that quenching star-formation 
(in the whole galaxy) is a slower process that causes also 
minor changes in their appearances due to disk-fading. 
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